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CN \ Abstract 

We consider cosmological dynamics in fourth order gravity with both f(R) and 
correction to the Einstein gravity {Q is the Gauss-Bonnet term). The particular 
case for which both terms are equally important on power-law solutions is described. 
These solutions and their stability are studied using the dynamical system approach. 
We also discuss the condition of existence and stability of a de Sitter solution in a more 
q-( general situation of power-law / and <3?. 

1 Introduction 

(TV 

Theories of modified gravity become very popular in the beginning of our century in attempts 
^ i to explain the accelerated expansion of our Universe which is proved observationally [TJ |2] . 

Another motivation is incorporating quantum corrections to General relativity (GR) which 
are believed to be present in the form of terms of higher order in the curvature. 

It is known that introducing additional terms, proportional to the Kretchmann invariant 
R^ uap R pV(Tp into the action often leads to solution with increasing spatial anisotropy [31 H] 
incompatible with the known picture of our Universe. This is the reason why modified gravity 
theories often deal with only /(i?)-corrections to Hilbert-Einstein action. However, special 
^ ; properties of the Gauss-Bonnet combination Q = R ta/as> R iaiap — AR^R^ + R 2 make f(Q) 

correction terms reasonable to study. Since Q is a topological invariant in 4 dimensions, terms 
linear in Q do not modify the equations of motion (if Q enters into the action in a product 
with function of another dynamical variable, like in string-gravity corrections, it modifies 
equations of motion and gives rise to interesting effects [SI El C])- Cosmological effects of 
f(R) correction terms have been studied recently for a number of different functions /. 

In the present paper we consider both f(R) and corrections to the Einstein gravity, 
find a situation where they are equally important and outline differences of cosmological 
behaviour in such combined model from models with only f(R) or $((?) terms in the action. 

Viable cosmological models have been constructed for each theory separately (see, e.g., [HI 
EJQIJJQTJ for Furthermore, finite-time singularities ([12]), cosmological perturbations 

( [13J ) , stability around a spherically symmetric static space-time ( [14J ) , and the number 
of degrees of freedom and ghost problems ([E]) have already been investigated. Recently 
more general theories with JF[R,Q) have become a subject of investigation [T4} [JB|. The 
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review [T7] considers some more complicated versions of modified gravity including non-local 
theories and solutions which can be of interest from the viewpoint of cosmic acceleration have 
been described. Our goal is different: we try to consider the general picture of cosmo logical 
dynamics in modified gravity theories applying the Dinamical system approach, used actively 
in f(R) gravity [HI [191 EO]- We make the next step taking into account the Gauss-Bonnet 
term restricting ourself in the present paper to a particular case where the function J-{R, Q) 
is a sum of a function of R and a function of Q . We find asymptotic regimes and studying 
their stability in such a modified Gauss-Bonnet gravity theory. 



2 Basic equations 

The action we study is given by0 

f(R) + HG) 



16nG 



N 



-gd 4 x + S M , (1) 



where f(R) and $(G) are general differentiable functions of the Ricci scalar R and the Gauss- 
Bonnet invariant G, Gn is the Newton constant. Variation of this action with respect to g^ v 
leads to the following field equations: 

Rnu — T^fiV — = SvtGatT^ , (2) 

where 

=V M V,F - g^DF - 2 9fiV RD^ + 2i?V M V^ - 4f2*V A V„f - 

- 4i#V A V M £ + AR^DZ + Ag^R^VaVpt + AR^VV^ - (3) 

- \g»v{FR + f(R) - $(£)) + (1 - F){R^ - \g^R) . 
The stress-energy tensor is 

Note that here and henceforth we use the notations F(R) = df(R)/dR, £((?) = d${G)/dG- 
As the attempt to consider the general f(R) even without the Gauss-Bonnet term [18] 
appears to be not fully successful and was criticized in [20] , we consider a particular simplest 
and mostly natural family of power-law functions in the Lagrangian. Namely, we study the 
case of f(R) = R + otR n and <&((?)= (3G m because it is important to reproduce GR results 
in a low-curvature limit. 

In the high-curvature limit, the Einstein contribution can be neglected, so we can choose 
f(R) = R n . The Gauss-Bonnet term is a total derivative in (3 + 1) space-time, therefore, a 
linear function of G in the action will not contribute to the field equations. We consider the 
Friedman-Lemaitre-Robertson- Walker metric for a flat isotropic universe: 

ds 2 = —dt 2 + a 2 {t)dx i dxi ; (5) 

using the Hubble parameter H = a/a, we rewrite the field equations in the following form: 

3H 2 = h\{FR + £Q - f - $) - 3HF - 12H 3 £ + 8nG N p] , (6a) 
b 2 

F = 4£H 3 - 4£F 2 + FH — 8£HH - 8nG N (p + p) - 2FH . (6b) 



1 The signature of the metric is assumed to be (— , +, +, +) and the speed of light c is taken to be equal 
to 1. 
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Moreover, we will need the continuity equation, the cosmological equation of state and 

expressions for the Ricci scalar and the Gauss-Bonnet invariant in terms of the Hubble 
parameter: 

= p + 3H(p + p) (7) 

R = Q(2H 2 + H) (8) 

g = 2AH 2 (H 2 + H) (9) 

p = up (10) 

3 The existence of de Sitter solutions 

The existence of de Sitter solutions, corresponding to an exponentially increasing scale factor 
(H — Hq — const), is an important opportunity in modified gravity theories. 
Taking the trace of Eq.(T2]), and assuming F = const, £ = const, T = we find: 

F(R)R + 2£(g)g-2f(R)-2$(g) = 0. (11) 

Consider the cases, mentioned in the previous section: 

1 f(R) = R n ,$(g) = f3g m : 

(n - 2)R n + 2/3(m - l)g m = . (12) 
Using the relations R = R = 12 H$, g = g = R^/6, we get: 

(n-2)i$ + 2/3(m-l)f-^J =0. (13) 

We can see that de Sitter is absent for the following combinations of m and n: 

1. n = 2 for any m 

2. m = 1 for any n 

3. n^2 ,m = n/2 

Thus de Sitter solutions exist for all other combinations of m and n. 

2 f(R) =R + aR n ,$>(g) =(3g m : 

2/3(m - l)g m + a(n - 2)R n - R = (14) 

and 

R 



2p(m-l)l-£J +a(n-2)R%-R = 0. (15) 

We can see that de Sitter is absent for the following combinations of m and n: 

1. n = 2, m = 1 

2. m = 1/2 and n = 1,2 

Thus, de Sitter solutions exist for all other combinations of m and n. 

We can see that adding $((?) function into the modified gravity action allows us to obtain 
new de Sitter solutions which do not exist in f(R) and pure Gauss-Bonnet gravity. 
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4 Cosmological dynamics 



Since the field equations in Modified Gravity theories are non-linear fourth order differential 
equations, which are very difficult to solve, we use the Dynamical system approach. It 
provides a powerful and relativity simple scheme for obtaining asymptotic solutions and 
investigating their stability 

The main goal of this approach is to get some autonomous system of first order differential 
equations, as a consequence of the cosmological equations of motion. In this approach, the 
dynamics of the Universe corresponds to motion along a phase curve, and stationary points 
represent some asymptotic regimes of the Universe evolution. 

Dividing Eq. ( |6a|) by 3H 2 , we find: 

., = J_ JQ_ _ J_ 5 F 4Hj 8nG N p 

6H 2 6FH 2 6FH 2 6FH 2 FH F 3FH 2 ' 1 ' 

Now we introduce the following normalized expansion variables: 

X =FH> ( 17a ) 

(17b) 



QH 2 
f 



6FH 2 



(17c) 



n = ^, d7d) 
- iG (17e) 



T 6FH 2 ' 

K = WW" < 17f > 

<l> = ^ ■ (17g) 
We take derivative of these variables with respect to the dimensionless time N = | lna(t)|: 

dx F 2 H 

dN = FH 2 ~ X - X JT 2 > (18&) 

d y R o h haw 

dN = 6lp- 2y H- 2 > _ (18b) 
dz R H 

IN = 6IP- ZX - 2Z IP' (18C) 
dQ 8nG N p H 

dN = WJp- nX - 2n H- 2 > (18d) 

^ = y + ^-^, (18e) 

dd) g g h , x 

W = *2i^ + *0H-*'- 3 *W (18f) 

^L = ^- kx -2k§- 2 . (18g) 

Using the new variables, Eq. (jl6p can be rewritten as 

\ = yJ r( f ) - Z - K - X -^J r Sl. (19) 
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We have obtained a set of seven equations, but our theory is of the fourth order. As we 
will see further, our dimensionless variables are not independent, consequently, this system 
is overdetermined, and our next purpose is to reduce its dimensionality 

To complete transformation, we need to express all right-hand side terms in Eqs. (|18[) 
through the variables defined in Eqs. (17). 
Using Eqs. ©,©,©, and (HD), we find: 



H 

R xy 



T^ = 2/-2, (20) 



6H 3 b ' 
87iGnP 
3FH 3 

G 



(21) 

3fi(l+w), (22) 



QH y 



y-1, (23) 

(24) 



24# 4 

Q (2(y-2) 2 + f) 



where 



, _ IrrR f _ d 2 f(R) , . 

If the functions f(R) and $({?) are chosen, we can close the system due to additional relations 
between our variables. 

The case of $((?) = /3G m . Using Eqs. © and Q, we get: 

= mK , (26) 
0(m-l)(2(y-2) 2 + f ) , s 

^ = - — ■ (27) 

Thus, using these relations with the constraint ( EES]) , we exclude n,ip and x from our au- 
tonomous system. As a result, we obtain a set of four equations: 

(28a) 

£L = 2L-zx-2z(y-2), (28b) 



In 



30(1 + u) - ttx - 2Q(y - 2) , (28c) 
<px - 2<p(y - 2) , (28d) 



m(2(y-2) 2 + f) 



rfiV (y - 1) 



where 



n-l + y- Z + ( p( m -l)(±--^- 



2(y-2) 



2 , 



x = ■—, , (29) 

i I #(m-l) v ' 

1 "I" &(j ,_i)2 

6 = IrrR _ (30) 
We would get additional relations between the variables for every particular function f(R). 
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5 The high- curvature case of f(R) = R n 

First of all, we can easily find in this special case 

y = nz, (31) 
b = n — 1 = const . (32) 

We can substitute these relations into f l28|) and exclude z from our system: 

* -2L-2,,fe-2), 
dJS n — 1 

^ = -3ft(l + w) - fiar - 2fi(y - 2) , 

# , m(2(y-2) 2 + ^ I ) (33) 
^ = — x ^-20(y-2), 

n_l + y 2^ + 0( TO -l)(±-?&=^ 



1 + 



y<f>(rn-l) 



(n-l)(y-l)2 

Note that the factor n — 1 is a denominator in our system. It means that the selected set of 
variables is useless for Gauss-Bonnet gravity with the Einstein term, where n = 1. 

Solving the system with vanishing left-hand sides, we find the following fixed points 
{x,y,Vl,<f>): 

(-1,0,0,0), 
(1 - 3w,0,2- 3w,0) , 

i g (n-l)(l+cj) 4n-3(l+m) 8n 2 -13n+3+qj(6n-3)(n-l) q\ 
V d n ' 2n ' 2n 2 ' U J ' 

/ 2(n-2) (4n-5)n n ^ 

\ 2n-l ' (ra-l)(2n-l) ' u ' u 2 ' 

(0,2,0,-§^), 
(-8m + 4, 0,0, m(8m ~ 5) ^ 



' "> (m-l)(8m-l) / ' 

The first four points are precisely the same as those obtained by Carloni, Dunsby, Capozziello 
and Troisi for i? n -gravity [19] Jj The Gauss-Bonnet term does not contribute here. Except 
the i? n -gravity solutions, here are two new stationary points: 



1. (0, 2, 0, It is the de Sitter point. 

2. (—8m + 4, 0, 0, (^z^^iiy )- The scale factor is a(t) = ao(t — to) 2 ", so this point repre- 
sents the well-known scale factor evolution dynamics of /(-R)-gravity, though the set 
of coordinates for this point is different from the first and second points, which also 
correspond to a ~ t*. 

It was shown in [21] (where combined effects of f(R) and RDR corrections were studied) 
that using the dynamical system approach we can obtain some new non-trivial solutions 
mainly where the correction terms in the action are equally important for a power-law 
evolution of the scale factor. Also there exist some additional relations between variables in 
this case. As we see further, the similar picture exists in our system. 

The case m = n/2 in R n + (3Q m models. There is an additional relation between <p 
and y in this special case: 

*=%V^- (34) 



1 Note that they have used a bit different set of variables and signature from our paper. 
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Consequently, one can reduce the system to two equations: 



dy 
dN 



x 



n-1 | l(f- 1)2^(^-1)* 



S-2 



dN 



y ?>2y n - 1 
-fi(3(l+u;) + 2(y-2)) - 



- 2y(y - 2) , 



(35) 



O 14-, y 2f (y-l)f n 2 2(y - 2) 2 



x 



1 + 



/3(f-l)2f^(y-l)f- 

3%yn-2( n _ 1) 



2 

-l 



It is useful to describe our future strategy. We have started from the model R + aR n + /3Q n 
and consider its high-curvature limit, where R can be neglected. We have two coupling con- 
stants, a and j3, however, when neglecting the Einstein term in the action, effectively our 
result depends only on their ratio j3/a. Thus we can assume a — 1. Next, we select the 
model by fixing n and find stationary points for some values of (3. Having this information 
we will discuss the situation in the general case. 

We mostly investigate the case of positive integers n and m, so n is chosen to be even. 
If n — 2, m — 1 so the Gauss-Bonnet term does not contribute. Thus, we study the case of 
n = 4 in detail. 

The case of n = 4,(3 = 1. Stationary points in this model, which have been found 
numerically are, by the scheme (x, y, Q, 0): 



CI 



9.358, 0.440, 0, 0.816), 
m:(0.576, 2.096, 0, 0.028), 



9. , 13 



3 / , 



1 15309^ 4 - 179820a; 3 +628974^ - 168876u- 1599187 2(0.625-0. 375u) 2 
288 (3w-13) 3 ~~' " > 9(1.625-0.375w) 3 



The corresponding eigenvalues are (for the last point eigenvalues depends on equation-of- 
state parameter u, so the values are presented for several particular co): 



Al : (7.918, 
Bl : (-3.672, 

f 

co = —1, 
co = -0.5, 
C\ : { oo = 0, 
co = 0.5, 

CO = 1, 



9.478 - 3co) , 
3.768 - 3co) , 
0.722, 
1.965, 
3.081, 
4.115, 
5.083, 



- 3.722 

- 3.653 

- 3.456 
-3.177 

- 2.833 



These points are summarized in the Table 1. 

The case of n — 4, /3 — — 1. Stationary points in this model are, by the scheme 
(x,y,Q,(j)): 



A2 :(-15.754, -0.625, 0, 2.396), 
B2 : (0.566, 2.094, 0, -0.028), 



C2 : 



9, , 13 



1 15309cj 4 - 179820^ 3 +596718u 2 - 159660^-1524947 2(0.625-0.375a>) : 



288 



(3w-13) 3 



9(1.625-0.375w) 3 



)■ 
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Table 1: Stationary points for R 4 + Q 2 



Point 


Coordinates (x, y, ft, 0) 
of stationary point 


Stability 

type 


The scale 
factor, a{t) 


AX 


x = -9.358 
y = 0.440 

ft = 
= 0.816 


Repulsive node 




Bl 


x = 576 
y = 2.096 

ft = 
= 0.028 


Attractive node 


Oo(t — to)~OA96 


CI 


9 9 

y 8 8^ 

o _ 15309l*; 4 -179820w 3 +628974w 2 - 168876oj - 1599187 


Saddle 


a (t-t ) 3(w+1) 


288(3u-13) 3 

/ _ 2(0.625-0.375w) 2 
^ ~~ 9(1.625-0.375w) 3 



The stability types for all this points did not change as compared with the case (3 = 1. The 
situation in all these cases is the same as has been found in R n gravity - only one stationary 
point (which corresponds to a phantom solution) is stable, others are unstable. 

The general case of n = 4. Now we try to describe what happens for an arbitrary 0. 
Substituting a(t) = a (t — t )A into Eq. fl6b|) . we get for R + aR 4 + (3Q 2 {n = 4) an algebraic 
instead of differential equation: 

' : (248832aA - 414720aA 2 + 214272aA 3 - 36288aA 4 + 



AT(t-t y 

+ 768/3 + 4608/M - 5376/3 A 2 + 27648a - 1 J 

- 2A\t - t ) 6 + 87cG N pA 7 {t - t )\l +u)) = 0. 

Neglecting matter and the Einstein term, assuming a = 1, we find the equation for A as a 
function of (3: 

- 144 - 1296A + 2160A 2 - 1116A 3 + 189A 4 - 4/3 - 24/3A + 28/3A 2 = . (37) 



A Stability analysis provided for a set of particular /3 reveals the following picture: 

• A<0. In this area there are only phantom branches, which are always stable. For 
/3 = ±1 they are solutions Bl, 2. Note that in the coupling constant interval —36 > 
/3 > —78.944 we do not find any phantom solution at all. 

• A>1. Branches of this region are always unstable. 

• AG [0, 3/7 w 0.428]. For this interval (3 e[-36,» -42.402] the branch stability depends 
on the parameter u (e.g., for (3 = —36 the eigenvalues are [—3, —3 — 3a;]). Considering 
also the matter-dominated solution, we get the following picture: there is some critical 
value 

8A-3 

u* = . (38) 
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Figure 1: Power index A of a vacuum power- law solutions depending on coupling constant 
for C = R 4 + (3Q 2 . Stable branches are solid bold lines, branches with stability depending 
on uj are dashed bold lines. 



In the case of uj < go*, the matter-dominated solution is stable while the vacuum power- 
law solution is unstable. In the opposite case uj > uj*, the solution with a ~ 

t l,A 

IS 

stable while a ~ £ 8 /[ 3 ( 1+w )] is unstable, i.e. solution with a larger power index is always 
stable. Thus, there is always one stable non-phantom solution for any uj in the range 
of the coupling constants -42.018 < (3 < -36. In the range -42.402 < (3 < -42.018 
(between the local maximum and minimum of function (3(A), see Fig. 2) we can find 
two stable vacuum solutions for some values of uj. 

A> 3/7 ~ 0.428. Here the vacuum power-law solution becomes unstable for any uj, 
while the matter-dominated solution a ~ t 8 /I 3 ( 1 + w )] i s still stable for some uj values 
close to —1. Furthermore, it becomes also unstable for (3 < —54, A > 0.791. 

However, it does not mean that there are no stable solutions for this range of coupling 
constants. In this section we considered models without the Einstein term, i.e., high- 
curvature regimes. Actually, earlier we noted that we should add R contribution 
for obtaining a late-time standard cosmology. This contribution gives us a de Sitter 
solution, which is stable in this range of coupling constants (see below). 



The case of general n. Note, that Eq.( l36p can be simplified only when the power-law 
ansatz terms originating from / and $ contribution are equally important. In this section 
we considered the case n = 4 for R n + [3Q% and discussed it in detail. Our studies show that 
the picture obtained does not change qualitatively for any n > 2. In all these models we can 
consider f3 as a function of the power index in the time dependence of the scale factor A, 



0{A) = -2 



(A - 1)(-24(A - !))-§ (-6{A - 2)) n (n - 3nA - 2 + A + 2n 2 A) 



(n - 5nA -2 + 2A + 2An 2 )(A - 2) 



(39) 



and admit some common features: 



1. There is a stable non-phantom solution in the coupling constants range (3 m i n < (3 < (3(0), 
where f3 m %n is the value of the function (3(A) at the local minimum point 



A 



1 — n + \/YJn 
2~ 



24n + 8 



ri 2 — 3n + 1 
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ca -44- 



-8.1 0.1 0^2 o!3 OA 0^5 0.6 0.7 

A, Power index of scale factor 

Figure 2: One of the branches plotted in Fig.l is shown here in more detail . The notation 
are the same as in Fig.l. 



P(A = 0) = —(6) a. Moreover, there is some critical value ui*, determining stability in 
this range: 

2nA-3 , . 

u* = — ^— (40) 

Solution with a larger power index (it is a matter dominated point or a vacuum power- 
law solution, depending on to) is stable. 

2. After A reaches the value A = 3/(2n — 1), this branch becomes unstable. 

3. There also exists a range of the coupling constants 



^-2 n^-Sn+l > < * < W)' 



where we have no phantom solutions at all. 

The only qualitative difference for the case of odd and non-integer n is the absence of .A-type 
solution in the range A > 1. 

We should conclude that for = —(6)2 there is a stable solution with A = 0. It means 
that this solution cannot be of power-law type. Using Eq. ( II 3p we easily find that it cannot 
be a de Sitter point. Since in Quantum field theory coupling constants are usually running, 
we do not consider in detail this solution existing only for a particular value of the coupling 
constant. 

It is interesting that in i? n -gravity with n > 2 we always have a stable phantom solution. 
By using the power-law ansatz it is simple to obtain that in pure Gauss-Bonnet gravity with 
J-"(R, Q) = Q m there is a phantom solution when m > 1/4. So, a combined effect of both R n 
and Q m correction terms may lead to solutions that are qualitatively different from those in 
pure R n or Q m theories. 
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6 Dynamics in R + (3Q m theories 

Now we consider cosmological dynamics in a theory where the only correction to Einstein 
gravity is a function of the Gauss-Bonnet term. It was shown in Sec. 5 that the set of 
dimensionless variables (TIT]) does not allow us to study this case, so we choose another set 
of variables. Substituting f(R) = R into (Q, we find the field equations: 

3H 2 = h&G-Q) -l2H 3 £ + 8nG N p, (41) 

= 4£# 3 - 4£H 2 - 8(HH - 8tvG n ( P +p)-2H. (42) 
We rewrite these equations in a dimensionless form: 

, iS * iHi+ ^. (43) 



6H 2 6H 2 s 3H 2 
^-4|-8^4-?^M-2|. (44, 



Substituting $ = (3G m , & = m/3g m , we get: 



The new dimensionless variables are 

<f> = ^ > ( 46a ) 

V> = 4#£ , (46b) 
8irG N p 

n = -^r , (46c) 

* = ^ ■ (46d) 

Taking derivative of the introduced variables with respect to N = |lna(t)|, we get the 
following set of equations: 



dQ 

dN 

dip 

dN 



-3fi(l + u) - 2Q(T , (47a) 
1ht + 4£, (47b) 



dN H 3 
Taking into account the relations 



m= m *$5- 2 * ir > (47c) 

d(J * 2a 2 . (47d) 



= ^-4£-2^a-3ft(l + w)-2a, (48) 

1 = <p - ip + tt , (49) 

g ii>(i + o) 

(50) 



i2 



iJ 3 m0 



- 4a - 2a 2 , (51) 
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we find the resulting autonomous system: 



^ = -3ft(l +w) -2£la, 



_^ = (0_l + fi )(l + CT )_20a, 

da_ " ' ^ 2 

dN 



7710 



(52a) 
(52b) 
(52c) 



Equating the left-hand side of our system to zero, we find the stationary point (0 = = 
0, n — 0, a — 0), which corresponds to a de Sitter solution. 

As in the previous case, nontrivial results have been obtained for a particular value of 
the power index m. 



The case m — \. There is an additional relation in the case of m — \: 



4> = -P 



(7+1 



We can substitute this relation into our system and find 
^ = -3Q(1 +w) -2tta, 

dN 8 v 7 



(53) 

(54a) 
(54b) 



Using the technique applied in the previous section, we find stationary points for particular 
values of 8. For 8 — 1 they are, by the scheme (0,-0, fi, a): 

A3: (0,-1, 0,-1), 
B3 : (-1.263,-2.263, 0, 8.582), 



. (_1 _1 _ , , V-9-3cj(1+cj) 1 6(3cu+l)-(5+3cu)y / -3-9^ _3 _ 3, ,\ 
L " 3 ' I 2 V 3 W ' 3w+l '6 3w+l ' 2 2 W i' 



The corresponding eigenvalues are: 



^3 : (4, -3u - 1) , 

B3 : (-11.582, -3u- 20.165), 

w = -l, 2.505, - 5.505 ^ 
w = -0.8, 1.996, -4.696 
w = -0.6, 1.345, - 3.745 
^ w = -0.4, 0.423, - 2.523 J 



C3 : 



> . 



Now turn to the case of 8 — — 1. We get the fixed points by the scheme (0, ip, Q, a): 

A3: (0,-1, 0,-1), 
V 3 : (0.2 63, -0.736, 0, -0. 582) , 

co . ( V-3-9u V-9-3cj(1+cj) 1 6(3^+l)+(5+3^) v / -3-9cj _3 _ 3 \ 
C0 ■ v 6 ' 3w+l ' 6 3w+l ' 2 2 W / ' 
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with the eigenvalues 



1.5 + 1.883z ' 
1.350 - 0.525i 
2.433 
2.296 



We summarize our results in Table 2. 

Table 2: Fixed points for R + f3^g 



Point 


Coordinates (0, ip, fl, a) 
of stationary points 


Stability 

type 


The scale 
factor, a(t) 


R + VG, P = l 


A3 


= 

i> = -1 

n = o 

a= -1 


Repulsive node 


a>o(t - 1 ) 


B3 


= -1.263 
V> = -2.263 

n = 
<j = 8.582 


Attractive node 




C3 


- 

1, _ V-9-3w(l+w) 
^ 3w+l 
q l 6(3u>+l)-(5+3w)V-3-9w 
" — 6 3w+l 

,-r 3 3,, 


Saddle 


ao(t- to) 3 < 1+w) 


R-VG, = -l 


A3 


= 
^ = -1 

n = 

(7= -1 


Repulsive node 


a (t - 1 ) 


V3 


= 0.263 
^ = -0.736 

n = o 

a = -0.582 


u < -0.611 
Saddle 
-1/3 > u > -0.611 
Attractive node 


ao(t — to)o^ 


S3 


- ^- 3 - 9w 

V 3w+l 

q 1 6(3u>+l)+(5+3u>)V-3-9w 

" _ 6 3w+l 

,-r 3 3,, 
° = ~2 ~ 2 U 


u < -0.777, Attractive focus 
-0.777 < u < -0.611, Attractive node 
-0.611 < u < -1/3, Saddle 


a (t-i ) 3(1+w) 



^3 :(4, -3u> - 1) , 
D3:(-2.417, -3a; - 1.834) , 



S3 : 





-1, 


—1.5 — 1.883z, 


1" 


-0.8, 


-1.350 + 0.525i, 


LO = 


-0.6, 


0.033, 


U = 


-0.4, 


0.196, 
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Using the power-law ansatz a ~ (t — t )A with the field equations f l4"Tj) we can investigate 
the general case: 

1 



: (6\/l - A + (3VQ(A + 1)) = . 



(55) 



2A 2 VT zr A(t - to) 2 

In this way one can express the coupling constant (3 as a function of A (see Fig 3.) Incidentally, 




-3 _2-1 1 

A, Power index of scale factor 



Figure 3: Power index A of a vacuum power- law solutions depending on coupling constant 
for R + f3\/Q. The notations are the same as in Fig.l 



when (3 > or f3 < — \/6, only phantom solutions are stable. When —a/6 < (3 < 0, there is 
only one stable non-phantom solution. It is either a vacuum power-law solution or a matter 
dominated point. 

Similarly to the situation considered in the previous section, there is a critical value 

2A -3 

u* = — ^— (56) 

that determines the power-law solutions' stability in this coupling constant range. It means 
that in this theory, like that considered in the previous section, there is a coupling constant 
range in which we have no vacuum phantom solutions. 

It should be noted that to consider a matter source with w > —1/3 we have to change 
the sign of the Gauss-Bonnet term, i.e., to consider a theory with R + [3y/—G- In such a 
theory we have a stable matter dominated point corresponding to a ~ £ 2 /[ 3 ( 1+w )l. 



7 Stability of de Sitter solution 

Let us discuss the stability of de Sitter solution. It is the case in which one can easily obtain 
a stability condition. Consider linear perturbations about the de Sitter point: 



H = H + 6H, (57) 

R = R + 6 (AH 5H + 5H), R = , (58) 

g = Go + 24(4H*5H + Hl5H) , Q = 24# 4 , (59) 

5R = 6(AH 5H + 8H) , (60) 

5G = 24(4H 3 5H + H$6H) . (61) 
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Substituting these relations into Eq. f lSaj) and keeping only linear terms with respect to 5H, 
we find: 

5H + 3H 5H + HUH ( - R — - 4 ) = (62) 

We are seeking solutions in the form 5H = Cie x+t + C2e x ~ t , and substituting it into Eq. (l62|) . 
we obtain an algebraic equation with the following roots: 



A± " ± f - \ {4&H^H^+mfMH ) " 4 1 1 (63) 

This shows that the de Sitter point is stable under the condition: 

48ff 8 $ 6< ;(ff<,) + 3g 2 /,R fl (g ) 1 

MW, < 4 (64) 

in agreement with [IB]. For f(R) = R + aR n , $((?) = we easily get: 

g2w -2 /M" - 2)2t12t-i + q W (n - l)^"' 1 

^° 1 + anl2^H 2 n - 2 < 1 ' ^ 

Substituting H from the equation of motion (ITS]) , we get the following stability condition 
of the de Sitter solution: 

- < -(65), n>2, (66) 
a 

->-(6i), l<n<2, (67) 
a 

- < -(6?), n<l. (68) 
a 

Comparing this conditions with the stability situation for power-law solutions in R n +f3Q^ 
we can make a conclusion that there is some non-trivial connection between these solutions. 
When the de Sitter point becomes stable, stable power-law solution become non-phantom. 
When this power-law solution loses its stability, the de Sitter solution remains stable. 



8 Conclusions 

In this paper, we have considered the cosmological dynamics in modified gravity with both 
f(R) and the Gauss-Bonnet term in the Lagrangian. In this theory, the case of R n + j3@2 
appears to be special one. In this case the terms in the equations of motion, originating 
from these two different modifications of the action, are equally important for power-law 
solutions. The Dynamical system approach works well only in this case and helps us to find 
exact solutions and study their stability. It was obtained that in the presence of Gauss- 
Bonnet term the general picture differs from that found for f(R) and /(C?)-gravity. There 
is a range of coupling constants in which there are stable non-phantom solutions for any 
integer value of n > 2 as well as for n = 1. 

It has been shown that there are some difficulties when the dynamical system approach is 
applied to f(R) + $(£?)-gravity. For power- law f(R) and $({?) studied in this paper it works 
well only while both terms are equally important on power-law solutions for the scale factor. 
It is useless for transient regimes, where only one term dominates. Thus we cannot investigate 
the Universe dynamics using only this method because different terms can dominate at 
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different stages of its evolution. It should also be noted that some cosmological solution can 
be lost when transforming the initial system to a set of first-order differential equations in 
expansion normalized variables, and the choice of these variables may be sometimes tricky 
(in our paper we use two different system of variables to consider n = 1 and n > 1 cases). 

We have also studied in detail the existence and stability conditions of a de Sitter solution 
for general power-law functions f(R) and $(£). 
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